Astronomy & Astrophysics manuscript no. berdyugina 


©ESO 2011 


September 15, 2011 





Upsilon Andromedae b in polarized light: New constraints on the 

planet size, density and albedo 

S. V. Berdyugina 1 , A. V. Berdyugin 2 , and V. Piirola 2 

1 Kiepenheuer Institut fur Sonnenphysik, Schoneckstrasse 6, D-79104 Freiburg, Germany 
e-mail: sveta@kis . uni-freiburg . de 

2 Finnish Centre for Astronomy with ESO (FINCA), University of Turku, Vaisalantie 20, FIN-21500, Piikkio, Finland 
e-mail: andber@utu.fi, piirola@utu.fi 

Received accepted ... 

ABSTRACT 

Context. Polarimetry is a novel tool to detect and characterize exoplanets and their atmospheres. It provides unique constraints on 
the planet orbit parameters, such as its inclination, orientation in space and shape, as well as reflective properties and composition of 
the atmosphere. These combined with spectroscopic and photometric measurements can fully characterize the planet with respect to 
its size, density, temperature, composition, etc., even for non-transiting systems. 

Aims. Polarized scattered light from the non-transiting hot Jupiter v And b is measured to further constrain its orbit, mass, density, 
and geometrical albedo. 

Methods. We obtained polarimetric measurements in the UBV bands over the orbital period and deduce an average peak-to-peak 
amplitude of (49 ± 5) x 10~ 6 in both Stokes q and u. 

Results. From our data we evaluate the orbit inclination ; = 111° ± 11°, longitude of the ascending node CI = 236° ±12° (or equiv- 
alently 56°), the effective size of the scattering atmosphere in the optical blue of 1.36 ± 0.20 Rj. These combined with spectroscopic 
measurements result in the planet mass 0.74+0.07 Mj, mean density 0.36±0.08 g cirT 3 , and surface gravity ~ 10 3 cm s~ 2 , which favor a 
close similarity of v And b to other inflated hot Jupiters. We also significantly improved the periastron epoch T p = JD2, 450, 032.45 1 , 
interior conjunction epoch T t = JD2, 450, 034.668, and periastron longitude a> = 279° ± 14°. The latter indicates that the apsidal 
resonance known for planets c and d includes also planet b. Obtained limits on the wavelength dependent geometrical albedo (average 
0.35) indicate its similarity to Neptune with peak reflectivity in the blue. Combinig all available measurements at various passbands, 
we construct a unified wavelength dependent albedo of an average hot Jupiter. It appears to be largely shaped by Rayleigh scattering 
in the blue and atomic and molecular absorption in the optical and near infrared. 

Conclusions. Our findings demonstrate the power of polarimetry for studying non- transiting exoplanets. 

Key words, planets and satellites: detection — planets and satellites: atmospheres — planets and satellites: fundamental parameters 
— polarization — techniques: polarimetric - planets and satellites: individual: Upsilon Andromedae b 



1. Introduction 

Upsilon Andromedae b (v And b) is a short period (4.6 days), 
non-transiting hot Jupiter orbiting an F8 dwarf together with 
other three planets ( |Butler et al.1[TW7l[T999)|Curiel et al.||201 1) . 



Planets in such systems are often found close to resonance orbits 
which is highly interesting for our understanding of planet for- 
mation and evolution. For instance, in the u And system planets 
c and d are in apsidal alignment associated with a secular reso- 
nance (Chiang et al. 2001 1, while planet e is close to an external 
3:1 resonance with planet c ( Curiel et al.| 201 1] >. Astrometric 



measurements of planets c and d with Hubble Space Telescope 
combined with radial velocity variations revealed their very low 
orbit inclinations and high masses (McArthur et al. 2010 1. Their 
mutual orbit inclination of 30° appears to be at the limit for this 
system to be stable (e.g., Michtchenko et al. 2006 ; Barnes et al. 
[2TjTT] i. 

Photons interacted with or emitted by an exoplanet are a 
key source of information on its physical properties. So far, 
our knowledge relies mainly on transition spectroscopy and sec- 
ondary eclipse techniques and, therefore, is largely limited to 
transiting planets. Exceptions are v And b and a few other hot 
Jupiters, for which a variable infrared flux was measured di- 



rectly (e.g., Harring toret^[2Tj06][Crossfield et aL][20T0] i, and 
HD 189733b, directly studied with the help of optical polarime- 
try ( jBerdyuginaetaL] |2008] [20TT] hereafter B08 and Bl 1). The 
latter studies further constrained orbital parameters and revealed 
a strong wavelength dependence of reflecting properties (geo- 
metrical albedo) of the planet with a maximum in the blue. 
This was found to be strikingly similar to Neptune which ap- 
pears bright in the blue due to Rayleigh scattering and absorp- 
tion in the red (jKarkoschka} |1994| |Sromovsky"l |2005|l. Thus, 



polarimetric studies are essential for non-transiting planets as 
they provide an inclination of the orbit (hence, true mass), ra- 
dius of the planet (hence, density), and direct access to plane- 
tary atmospheres (Seager et aT] |2000| [Stam et alT 2004| Fluri & 



Berdyugina 2010 



Berdyugina} |201 1| . 



This paper reports measurements of the reflected light from 
a non-transiting extrasolar planet using polarimetry in the UBV 
bands. Form these data we evaluate the inclination and orienta- 
tion of the v And b orbit and the planet radius, which result in 
reliable estimates of its mass and mean density. In addition, we 
find indications of an atmosphere which efficiently scatters the 
stellar light in the blue (UB bands). Previous polarimetric obser- 
vations of v And b by Hough et al. (2006 ) were made in a broad 
red band (covering RI bands), and no detection or limits were re- 
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ported, which implied a very low albedo in the red. We estimate 
an average optical albedo of u And b of 0.35, similar to that 



of HD189733b (Bll ) and Kepler-7b ( Kipping & Bakosj |20TT 
|Demory et al.||20l"T| ). These three planets are inflated and highly 
reflective hot Jupiters, most probably due to Rayleigh scatter- 



ing on condensates in high altitude haze or clouds ( Berdyugina 
[20TT]|Demory et al.||20lT] l. 

2. Observations 

The observations were carried out in 2009 October 27 to 
November 3 and in 2010 September 15-22 with the TurPol 



broad-band (UBVRI) polarimeter (Piirola| [T973||1988| ) mounted 
in the Cassegrain focus at the 2.5m Nordic Optical Telescope 
(NOT). The observing procedure was the same as in Bll. In 
particular, the Stokes gw-system was defined in accordance with 
the common agreement that positive and negative q are in the 
north-south and east-west directions, respectively, while positive 
and negative u are at an angle of 45° counter-clockwise from the 
positive and negative q. 

To measure linear polarization (normalized Stokes qu), the 
polarimeter is equipped with the following elements along the 
optical path: (1) a half-wave plate retarder rotating at 22.5° in- 
tervals to modulate between polarization states, (2) calcite block 
to split the light onto the ordinary 'o' and extraordinary V paral- 
lel light beams, (3) two diaphragms to limit the light integration 
area around each beam, (4) chopper to quickly (25 Hz) switch 
between the ordinary and extraordinary beams, (5) dichroic mir- 
rors to separate the passbands, (6) field lenses to direct each of 
the two beams on the same detector, and (7) five photomulti- 
pliers as detectors for the five UBVRI bands. A drawing of the 
principle scheme for one channel can be found in Piirola ( 1973[ l. 

There are a few important advantages in this scheme which 
enable achieving very precise polarimetry with TurPol. 

First, since the calcite block is located above the diaphragms, 
an additional e-component from sky background, (sky) e , is su- 
perimposed to the (star+sky) fl component, while an additional 
(sky) a is sperimposed to (star+sky) c . Thus, both components of 
sky polarization pass both diaphragms. Since they are of the 
opposite polarization states, they cancel out exactly, indepen- 
dently on the amount and origin of background polarization. 
Therefore, polarization due to, e.g., scattering in the Earth at- 
mosphere, moonlight, or resolved extended sources (nebulae) 
is exactly eliminated at the time of target observations. This is 
a great advantage of TurPol in comparison with the polarime- 
ters PlanetPo l of |Hough"eTaLl ( |2006| and its copy POLISH of 
Wiktorowicz ( 2009} where a Wollaston prism is used instead of 
a simple calcite block. The prism splits the light with a large di- 
vergence angle between the ordinary and extraordinary beams, 
and therefore the sky beams do not overlap, and their polariza- 
tion does not cancel out. Therefore, such polarimeters should 
measure very carefully background fluxes in each beam sepa- 
rately to deduce the sky polarization and subsequently subtract 
it from target measurements. This is a significant source of errors 
in such schemes. For instance, |Lucas et aL] ( |2009] ) reported that 
dust blown from Sahara to Canary Islands can produce sky po- 
larization of order 10~ 5 at large zenith distances. The compensa- 
tion scheme of TurPol guaranties that this spurious polarization 
is completely eliminated at the time of measurements. 

Second, the rapidly rotating chopper (25 Hz) alternates in- 
tegration of the ordinary and extraordinary beams on the same 
detector. This strongly diminishes systematic effects due to in- 
strumental sensitivity (flat-field effects), seeing and variable sky 
transparency. Again, this is an advantage to other polarimeters 



where integration of orthogonal polarization states is carried out 
with different detectors and at much longer modulation times. 

Third, thanks to the dichroic mirrors and five detectors, sev- 
eral passbands can be measured simultaneously. This signifi- 
cantly increases the efficiency of observations and provides in- 
formation on spectral distribution of polarization, which is valu- 
able for identification of its source. 

Finally, each pair of the observed Stokes q and u was cal- 
culated from eight exposures at different orientations of the re- 
tarder, which allowed us to avoid systematic errors due to imper- 
fections of the retarder. 

To calibrate the data with high accuracy, we dedicated a sig- 
nificant amount of observing time for determining systematic 
effects. For this purpose, we observed 10 nearby, bright, non- 
variable, and nonpeculiar stars (V=5 m -6 m , d < 50 pc) of spec- 
tral classes F-G. This range allowed us to maximize the photon 
flux in the blue and observe without neutral density filters. The 
stars were observed on the same nights as v And at various par- 
allactic angles to determine instrumental polarization which in- 
cludes contributions from the telescope (P t ei) an d the polarime- 
ter itself (Pin). The former is rotating during observations due 
to the alt-azimuthal mount and the latter is constant. Total inte- 
gration time of 3-4 hours for each star allowed for accuracy of 
(l-2)xl0~ 5 . All stars, including scientific targets, were observed 
within ±30° zenith distance to minimize variations of instrumen- 
tal and sky polarization and to maximize the photon flux. In ad- 
dition, for calibration of the polarization angle zero point, we 
observed highly polarized standard stars HD25443, HD161056, 
and HD204827, with the accuracy better than 1°. 

The polarization calibration data are presented in Table [T 
f t ei corresponds to the telescope position angle zero. Rotationa 
transformation of these values are applied as a function of the 
angle of the telescope relative to the equatorial system, due to 
the alt-azimuthal mounting. 

After the re-aluminization of the primary mirror of the NOT 
in spring 2009 P te i reduced by about a factor of 10 as compared 
to 2008 (see Bll) and became comparable with measurement er- 
rors. In 2010 it slightly increased but still remained significantly 
lower than in 2008. This is also confirmed by measurements in 
June 2010 during another observing run with TurPol. P; n may 
include small effects from residual interstellar polarization of the 
standard stars, although these effects tend to average out for a 
sample of stars in different directions on the sky. In any case, 
uncertainties of P m will only cause a constant shift of the q- and 
w-curves, and have no effect on the amplitude of the variations 
we are using for our analysis. 

We observed v And by making 10 s exposures in the UBV 
passbands simultaneously, during two hours on each clear 
nights. The best photon statistics was in B because of the stel- 
lar spectral distribution. Since the star is very bright, we used 
a 10% neutral density filter which allowed the maximum pho- 
ton flux in B below the saturation limit. Because of bad weather 
we lost six nights. Errors of individual measurements obtained 
during one cycle of the retarder were conservatively assigned 
to either statistical or photon noise errors, whatever was larger. 
This may have lead to some overestimation of the errors, as the 
normalized x 2 of our fits to the data was about 0.8. By aver- 
aging individual data points we obtained nightly measurements 
(35 in total) with mean individual errors of 4.5 xl0~ 5 in U and 
V and 2.5xl0~ 5 in B. The data were corrected for the telescope 
and instrument polarization. The calibrated measurements are 
shown in Figure [T] The following peak-to-peak amplitudes in 
the individual UBV bands were evaluated (in the 10~ 5 scale): in 
Stokes q, respectively, 4.6±1.9, 5.8±1.0, 2.5±1.8; and in Stokes 
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Table 1. Polarization calibration data. Stokes parameters [q, u] ± <x are given in units of 10" 5 . 



Target 


Sp. 


U 




B 




V 




UBV 






October 2009 












HD24740 


F2IV 


[2.1,-2.7] ± 3.1 


[-1.3,-3.5] ± 


2.3 


[-5.7, 2.1] ± 


4.1 


[-1.0,-2.3] 


±2.2 


HD29645 


GOV 


[-5.0,-2.5] ±2.1 


[-0.7,-1.0] ± 


1.4 


[0.6, -2.8] ± 


3.2 


[-1.7, -1.6] 


± 1.6 


HD33167 


F5 V 


[-2.2, 3.2] ± 2.9 


[-4.2, 1.6] ± 


1.9 


[-1.8,-7.5] ± 


4.4 


[-3.4, 1.0] 


±2.2 


HD35296 


F8 V 


[3.9, 1.0] ±2.0 


[3.4, 1.8] ± 


1.4 


[2.8, 5.1] ± 


2.7 


[3.4, 2.1] 


± 1.4 






[2.5,-1.3] ± 1.2 


[-0.5,-2.8] ± 


0.8 


[6.0,-3.8] ± 


1.7 


[1.2,-2.5] 


±0.9 


P tel 




[-1.4,-2.3] ± 1.2 


[0.2,-1.1] ± 


0.8 


[-0.8,1.6] ± 


1.7 


[-0.4,-1.1] 


±0.9 


HD 16234 


F7 V 


September 2010 
[-4,-1] ± 3 [-2, 2] 


± 2 


[2, 0] 


± 3 


[-1.5,0.8] 


± 1.8 


HD29645 


GOV 


[9, -2] ± 5 


[-2, 1] 


±4 


[-4, 4] 


± 3 


[-1.0, 2.0] 


±2.8 


HD159332 


F6 V 


[-2, 0] ± 3 


[0, -5] 


± 3 


[-1,0] 


±4 


[-1.0,-2.0] 


±2.3 


HD 182807 


F7 V 


[-3, 6] ± 9 


[4, 14] ± 


10 


[-10,-4] i 


: 12 


[-2.3, 6.4] 


±7.3 


HD185395 


F4 V 


[8, 4] ± 5 


[6, 0] 


± 5 


[5,0] 


±6 


[6.5, 1.5] 


±3.8 


HD225239 


G2 V 


[1, 1]±4 


[1,0] 


± 3 


[3, -3] 


± 3 


[1.8,-1.0] 


±2.4 


Pin 




[1.5,-5.1] ± 1.6 


[-2.7,-2.9] ± 


1.4 


[0.0,-0.1] ± 


1.6 


[-0.6, -2.7] 


± 1.1 


P lei 




[-5.2,0.9] ± 1.6 


[-4.1,0.1] ± 


1.4 


[-0.7,0.1] ± 


1.6 


[-3.4, 0.3] 


± 1.1 



u 4.9+1.7, 6.1 + 1.1, 2.7+3.2. The uncertainties were calculated 
as standard deviations from the best fits (see Section|4]i. Thus, a 
statistically significant signal was measured only in the B band. 
However, by averaging altogether the f/BV-band data we can 
further improve the statistics and determine the Stokes q and u 
amplitudes, respectively, (48±4)xl(T 6 and (50+5)xl(T 6 (shown 
in the lower panels of Figure [TJ. This is the highest accuracy 
achieved so far in blue wavelengths. Even a conservative error 
of +7xl(T 6 (which is the mean error of the binned data) sug- 
gests the signal to be real. This average error is, however, not 
used in our subsequent interpretation and provided here only for 
evaluating the statistics of the data. A Monte Carlo analysis and 
parameter fits presented in Sect. |4]employ the individual nightly 
measurement errors. 

As an additional test for robustness of the detection, we made 
a model fit to the standard star data (Table [TJ phased with the 
v And b period. The maximum peak-to-peak amplitude found 
in both Stokes q and u was (2 ± l)xl0~ 5 , which is within the 
measurement errors and confirms that the evaluated amplitude 
of (5 ± 0.5)xl0~ 5 is not due to calibration errors. 



3. Possible sources of polarization 

We interpret the observed polarization as the result of scattering 
in the atmosphere of planet b which varies as the planet orbits 
the star. This is the innermost planet in the system, and scatter- 
ing polarization from the more distant planets in the system is 
expected to be significantly smaller than 1CT 6 . Two peaks per 
planet period in the phase curves suggest that the signal is due to 
reflection. Moreover, planet parameters determined from them 
agree with previous results (see Sect. [4]). 

There could also be a stellar contribution to the polarization. 
This can arise due to the Zeeman effect in spectral lines form- 
ing in magnetic stellar regions, such as starspots and plages. In 
order to estimate the amount of polarization due to active re- 
gions, we need both the stellar rotation period and the amplitude 
of variability. Measurements in the chromospheric Can H&K 
lines indicate a very low activity and variability of the star (e.g., 



Simpson et al. 2010 1. Because of that, its rotational period is 
poorly known. The Can H&K records reveal weakly significant 
per iods of 7.3 days ( Simpson et al.||2010"] l, 1 1-12 days ( Henry et 



al. 



al. 



2000; Shk olnik et al.||2008[ ), and even 18-19 days ( |Henry et 
2000} . For the period of 1 1-12 days to be due to rotation, the 
stellar radii should be about 2.1 - 23R a . This can be excluded 



-20 



-40 



10 
5 


-5 
-10 




iP 20 - 



-20 



V 



-40 ■ 



-0.4 



-0.2 0.0 0.2 
Phase 



0.4 



-0.4 



-0.2 0.0 
Phase 



0.2 0.4 





• 







10 r 

5 ; 

: 

-5 : 

-10 : 



-0.4 



-0.2 0.0 0.2 
Phase 



0.-1 



-0.4 



-0.2 0.0 
Phase 



0.2 0.4 



Fig. 1. Stokes q and u with ±l<x error bars for v And b. Upper 
panels: nightly UBV measurements for each band separately 
(squares). The U and V data are shifted in vertical by ±3xl0 -4 
for clarity. Lower panels: the UBV data averaged together. The 
mean error of the binned data is 7xl0~ 6 and the standard devia- 
tion is 4.5xl0~ 6 . Curves are best-fit solutions (Sect. ^ obtained 
by using orbital parameters listed in the last column of Table [2] 
for each band individually (upper panel) and all bands simulta- 
neously (lower panel). Phase 0.0 corresponds to the periastron 
epoch. 



based on the estimate of the stellar angular diameter of 1.1 mas 
from infrared flux measurements (Blac kwell et al.[[l990| ) and the 
parallax of 74 mas ( van Leeuwen||2007[ ), leading to the radius of 
1 .6R & , which is also in agreement with the value of 1 .64R sug- 
gested by|Takeda et al. (|2007[ l. A combination of vsku = 9.6 
( |Valenti & Fischer 2005) 1 and! .6R Q suggests the stellar rotation 



period to be shorter than 8.5 days. 

To obtain additional constraints on the period and to esti- 
mate typical projected area of active regions, we analyse here 
the Hipparcos photometry ( |Perryman et al. 1997K There 



93 measurements taken during 23 days distributed over about 
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Fig. 2. Hipparcos photometry of u And. (a) Original nightly av- 
erages (squares) and a double-harmonic fit (curve) revealing a 
possible long-term modulation of 684 days, (b) Periodogram of 
the data from panel a. (c) Nightly averages from which the long- 
term modulation was subtracted. The light curve is phased with 
the 7.88 d period, (d) Periodogram of the data from panel c. 



3 years (Fig. [2]). A double-harmonic fit to these data indicates 
a possible long-term modulation (684 days), while a Lomb- 
Scargle periodogram suggests a short-term variability of 7.88 
days (both periods are of low significance). Variations of the 
mean stellar brightness on time scales of several years are typi- 
cal for solar-type spotted stars (Berdyugina 2005). They can be 
caused by cyclic growth and decay of active regions as well as by 
their rearrangement on the stellar surface due to global magnetic 
field evolution. The short period is close to what is expected for 



the stellar rotation and to the period of 7.3 days by |Simpson et al.| 
(2010). It is better visible after removing the long-term modula- 
tion from the data (see lower panels in Fig.|2|i. If true, this period 
implies the stellar rotation axis inclination of 68°. Amplitudes of 
the long-term and rotational modulations are about 0.004 mag 
and 0.003 mag, respectively. We consider these as upper limits 
on stellar photometric variability. 

The photometric variability can be caused by either dark 
spots or bright plages (or both). We estimate that in either case 
the projected area of an active region causing such a variability 
is of the order of 0.5% (slightly depending on the temperature 
contrast). This is comparable on average with the solar variabil- 
ity. A linear polarization signal from such an active region in- 
tegrated over the stellar disk and a broad passband would be 
significantly smaller than 10~ 6 . Other effects, such as symme- 
try breaking of the stellar limb polarization due to active regions 
and stellar oblateness, are even less significant (see, e.g., esti- 
mates for HD189733 in Bll). Thus, we are safe to assume that 
the contribution from stellar magnetic regions is negligible. 



4. Orbital parameters 

Modeling observed variations of Stokes q and u enables re- 
construction of the orbit spatial orientation and an estimate 
of the effective size of the scattering atmosphere. Similar to 
the analysis by B08 and Bll, we employ a model based on 
the Rayleigh-Lambert (RL) approximation for an extended star 
(Fluri & Berdyugina 2010| hereafter, FB10), i.e. assuming 1) 
Rayleigh scattering for polarization and 2) the Lambert sphere 
with the geometrical albedo A g =2/3 for intensity. The second 



assumption has a negligible effect on the phase curve shape for 
unresolved systems (FB10). 

Since the albedo of the atmosphere is fixed in this model, 
the wavelength dependent properties of the planet are effectively 
included into the radius of the RL-atmosphere R RL , which rep- 
resents a geometrical limit for the unity optical thickness in the 
atmosphere in a given passband. This implies that a possible gas 
layer above the reflecting surface (e.g., clouds) has no influence 
on the flux or polarization (i.e., very optically thin). If the ra- 
dius of the planet is known from other measurements, e.g. from 
transits, the geometrical albedo can be evaluated as in B 1 1 . Even 
though these radii and albedos should only be considered as lim- 
its for a spherically symmetric case, they are useful for evaluat- 
ing fundamental properties of the planet. 

We used the^ 2 minimization procedure described and tested 
by FB10. Values of parameters previously evaluated from spec- 
troscopy were fixed. These were the orbital period P, periastron 
epoch T p , semi-major axis a, eccentricity e, and the radius of 
the star R t , which was considered to be a limb-darkened sphere 
2000). We considered three most accurate sets of pa- 



( Claret 



rameters obtained by Butler et al. (2006 ),|McArthuretal. 



(2010), 



arimet- 



(2C 

and |Curiel et al.| ( |20l'l) (see Table|2|i and searched for po 
ric solutions for each set separately. Our free parameters were 
the orbit inclination i, longitude of the ascending node Q, peri- 
astron longitude a>, and radius R^. Initially we also varied ec- 
centricity, but the amount of the data was only sufficient to con- 
clude that a non-zero eccentricity is necessary to achieve better 
fits. Therefore, we fixed the eccentricity to spectroscopic values. 
Since the orbit is nearly circular, the maxima are expected near 
the elongations. The eccentricity results in shifting the maxima 
closer to the periastron, which constrains a> and T p . The value 
of the latter was adjusted to place phase 0.0 at the periastron. 
The epoch of the interior conjunction (transit epoch T t in occult- 
ing systems) was assigned to the system brightness minimum. 
The deduced epochs agree well with those determined by others 
(Table [2). 

We emphasize that for non-transiting planets polarimetry 
provides a unique opportunity to directly evaluate the orbit in- 
clination i and Q. Moreover, it is possible to distinguish between 
inclinations smaller and larger than 90°, which is not possible 
from transit data. In our model (FB10) the inclination is de- 
fined in such a way that the planet revolves counterclockwise 
as projected on the sky for 0° < i < 90° and clockwise for 
90° < i < 180°. Further, Q. varies from 0° to 360° starting from 
the north and increases via east, south, and west. Note that there 
is ambiguity of 180° in the Q. value, which is due to the intrinsic 
property of Stokes q and u to remain unchanged under the rota- 
tion by 180°. The periastron longitude a> varies from 0° to 360° 
starting at the ascending node and following first the ascending 
part of the orbit. 

The four model parameters were first estimated from all 
UBV measurements simultaneously, since this provided the 
highest accuracy of the orbital solution (31 degrees of freedom). 
For each spectroscopic solution we found practically the same 
values of i, Q, and Rrl. The values of u> were found in agree- 
ment with the corresponding spectroscopic estimates. Their dif- 
ferences between the sets are due to differences in P, T t , and 
e. The solution with the smallest standard deviations was found 
with parameters by Mc Arthur et al. ( 2010[ l who combined spec- 
troscopic and astrometric measurements for the planets b, c, and 
d. Our solution is presented in the fourth column of Table|2] and 
the fits are shown in the lower panels of Figure [T] The errors of 
the parameters were estimated using Monte Carlo simulations 
as described by B08. The x 1 contours of this solution and re- 
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Table 2. Orbital and physical parameters of v And b. 



Parameter 


value 


value 


value 


value 


P, days 


4.01 11 1(6) 


4.61 / 1 1(1) 


4.61 /U3(3) 


A £\H 1 1 [51 

4.61 / 1 1 L ' 




1802.6(7) 






32.45 l(l) m 


1 t 


1802.97(3) 


34.1(3) 


5.37(5) 


34.668(1 ) m 


a, AH 


0.060(3) 


0.0594(3) 


0.0592217(2) 


0.0594 151 


e 


0.02(2) 


0.012(5) 


0.0215(7) 


0.012 [5) 


i° 








111(11) 


of 


63(7) 


44(26) 


325(4) 


279(14) 


£1° 








236(12) [8] 


M/Mj sin (' 


0.69(6) 


0.69(2) 


0.688(4) 




M/Mj 








0.74(7) 










1.36(20) 


p, g cirT 3 








0.36(8) 


g, 10 3 cm s -2 








0.99(46) 



Notes. 



[i] 



Butler et al. 



(20061 



[2] 



McArthur et al. 



1 2010 1; 131 



Curiel et al. 



J20TTI; 141 this paper; [5] fixed parameter; [6] JD 2,450,000+; m T p 



1800.805 and 4.749, T t - 1803.021 and 6.965, to compare with LI J and |3J, respectively; [8] £1 can also be 56° due to the 180° ambiguity. Errors 
corresponding to the last digits of the values are given in parentheses. 
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Fig. 3. The x 1 contours of the best-fit solution for v And b in the 
UBV bands simultaneously (solid lines). The three contours are 
shown for the confidence levels 68.3%, 90.0%, and 99.0%. Due 
to the intrinsic ambiguity in Q only one of the two equal min- 
ima separated by 180° is shown in each plot. The Monte Carlo 
sample solutions are shown with dots. 

suits of the Monte Carlo simulations are shown in Figure [3] The 
Monte Carlo tests were concentrated near the x 1 minimum for 
all parameters, which indicates that our fit to the data is robust to 
errors. 

The deduced epochs for planet b agree well with those 
by |Butler et al.] ( pU06] l and |McArthur etalj (pOTO) !, but dif- 
fer slightly from estimates by |Curiel et al. " f2011[ ) whose so- 
lution includes four planets and a trend but is based on spec- 
troscopy only. The inclination of the planetary orbit found from 
polarimetry is close to an indirect estimate of the inclination of 
the stellar rotation axis deduced from the stellar rotation periods 
and projected rotation velocity: 58°±8° from the period of 7.3 
days (Simpson et al.|[2~010[ ) and 68° from the period 7.88 days 
(Sect. [3J. In fact, their complementary values 180°-58°=122° 
and 180°-68° = 1 12° should be considered. Thus, the assumption 
that the stellar rotation axis is aligned with the normal to the 
planetary orbit made by Simpson et al. (2010 1 seems quite good 
for v And b, but is not true for the planets c and d in this system 
according to astrometric measurements (McArthur et al. 20 1 0) > . 
The found inclination also satisfies the lower limit i > 28° de- 



duced from infrared flux variations (Crossfi eld et al.| |2010) > and 
is outside the limits i - 90° ± 8° implying transits. It is interest- 
ing that Q=236°of planet b is the same as for planet c (237° + 7°, 
|McArthur et al.||2010) and that u of planets b, c, and d ( |Curiel et| 



JL 0- 



< 




2 0-2 
Aa (mas) 



Fig. 4. The orbit of v And b as projected on the sky. Solid and 
dashed lines indicate parts of the orbit in front of and behind the 
sky plane, respectively. The circle in the center depicts the star 
and the smaller circle on the orbit the planet. The direction of 
the orbital motion is indicated by the arrow. The reconstruction is 
made for the orbit parameters listed in the last column of Table|2] 
stellar radius R t /R = 1.631 (Baines et al. 2008), and distance 
d-YiAl pc. The planet is positioned at the periastron. If £2=56°, 
the parts of the orbit in front of and behind the sky plane are to be 
swapped together with the periastron position but the direction 
of the orbital motion is preserved. Positive directions of Stokes 
q and u and orientation on the sky are also shown. 



al. 2011 



are very close to each other too. This indicates that the 
apsidal resonance known for planets c and d includes also planet 
b, which provides further constraints on the dynamical history 
of this system. 

The full set of the orbital parameters of v And b allows us to 
depict its orbit as projected on the sky plane and to indicate the 
direction of the orbital motion. This is shown in Figure|4] Such a 
plot is useful for future direct studies of the system as the planet 
coordinates can be predicted quite accurately. 
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5. Physical properties of v And b 

Knowing the inclination of the orbit allows evaluation of a true 
mass of the planet 0.74Mj. This combined with the radius of 
1.36/?j leads to quite low mean density and surface gravity 
(Table [21, which identify v And b as an inflated, low density 



planet (e.g., Baraf fe~et al.| |2010 i. For instance, it is almost a 
twin to HAT- P-9b (|Shporer et al.| |2009| l and is very similar to 
HD209458b UMazeh et al.||2000), WASP- lb ([Collier Cameron| 
[etaL] |2007l >, and HAT-P-lb ( [Bakos et al.|[2007| ), to name a few. 
Interestingly, they all orbit F7-G0 dwarfs at 0.05-0.06 AU. Note 
also that about a quarter of known transiting planets have mean 
density lower than 0.4 g cirT 3 . Our 3<x upper limit on the surface 
gravity g < 2.3 x 10 3 cms -2 for v And b is consistent with the 
limit g < 2.1 X 10 3 cms -2 deduced from 24 fim flux variations 
( |CrossneldetaLl[20T0l i. 

These and earlier IR flux measurements by Harrington et al. 
(2006) revealed a large temperature contrast between day and 
night sides of the planet. Assuming R(24fim) = l.3Rj and zero 
albedo, Crossfield et al. obtained the contrast AT > 900 K and 
the lower limit on the orbit inclination of 28°. Using our inclina- 
tion value of 111°, we can evaluate the contrast AT ~ 1000K. 
The assumption on zero albedo at 24 jjm agrees with our conclu- 
sion that it diminishes strongly toward the red (see below). 

To investigate wavelength dependence of polarization and 
planet properties, we also modeled each UBV passband sepa- 
rately while having fixed all other parameters but the radius. We 
obtained R Rh = 1.34±0.35, 1.50+0.28, and 1.00+0.40 R 3 for U, 
B, and V, respectively. The corresponding fits are shown in the 
upper panels of Figure[T] Even though we obtained a statistically 
significant signal only in B, the other bands provide useful limits 
at the corresponding wavelengths. An estimate of the polarimet- 
ric signal in red wavelengths (RI bands) can be inferred from 
the upper limit by [Lucas et al. (2009 1 for t Boo b, when scaling 
it to the orbit size of v And b. This resulted in the polarization 
amplitude of < 3.3 x 10" 5 and Rrl(RI) < 0.36flj. 

Similar to HD189733b (Bl 1), the BVRI polarization ampli- 
tudes indicate the dominance of Rayleigh scattering in the op- 
tical, as they can be approximated with a Rayleigh law (Fig. [5] 
upper panel). The upper limit in U is however significantly lower 
than this relation and even smaller than in B. As pointed out by 
Bll, this can be a signature of additional opacity in the near 
UV which is perhaps stronger in v And b than in HD189733b. 
Taking into account that u And b is of lower density (more in- 
flated) and orbits a hotter and metal-rich star ([Fe/H] ~ 0.1), this 
is in accord with the suggestion by Burrows et al. (2007 2008 ) 



that inflated planets may have larger optical opacities (at higher 
altitudes) resulting in temperature inversion (stratosphere) and 
higher infrared fluxes. One possible source of such opacity is 
absorption in TiO and VO bands (e.g., |Fortney et al." 2008). Our 



data indicate that there may be an additional source in the UV, 
e.g., absorption by a haze of photolysis products or Raman scat- 
tering. Being confined to the UV, this additional opacity may 
play a significant role in heating of the planet. 

The large range of R RL in different passbands indicates a 
strong wavelength dependence of the geometrical albedo, as 
discussed by Bll. Assuming a maximum physical radius of 
the planet of 1.5 Rj, we can evaluate geometrical albedos A g = 
0.53+0.27, 0.67+0.24, 0.29+0.23, <0.04, in U,B,V,RI, respec- 
tively. These are plotted in Figure [5] (lower panel). Then, the 
average albedo for the range of 300-500 nm is 0.60+0.18, and 
for the entire optical range is ~0.3. Similar to HD189733b, we 
find a resemblance of A g (A) to that of Neptune, whose albedo 




1000 



Wavelength [nm] 



Fig. 5. Polarimetric amplitudes Aq (top) and geometrical albedo 
(bottom) measured in various passbands. Our new data are 



shown by open circles and the upper limit for r Boo b (Lucas et 
|aT)|2009 ) scaled to the orbit of u And b by a triangle. Horizontal 
bars show the passbands' widths. Upper panel: thin horizontal 
lines indicate Ag for different Rrl/Rj', solid (red) curve is the 
Rayleigh law Aq = 2.1 x 10 n ,T 4 scaled to fit the BVRI data. 
Lower panel: solid curve is the geometrical albedo of Neptune 



(Karkoschka Q994 ) smoothed by a lOOnm boxcar. Our data in- 
dicate the dominance of Rayleigh scattering in the planetary at- 
mosphere. 



is defined by scattering in the blue and by absorption in the red 
( |Sromovsky|[2005] ). 

A marginal spectral signature of the reflected light in the 
wavelength range 380-650 nm was reported by Collier Cameron 



et al. (2002) at a projected orbital velocity amplitude of the 
planet 132kms~ 1 . Assuming orbital inclination of 70°-80°and 



flat geometrical albedo of 0.42 (Class V models, Sudarsky et al. 



2000), they estimated the planet radius of 1. 34+0.1 7 R } . These 
assumed and deduced parameters appear very close to our re- 
sults. For instance, our orbital parameters predict the veloc- 
ity amplitude of the planet to be 138kms~'. Thus, the Class 
V models with highly-reflecting clouds at lower pressure alti- 
tudes apparently hold for the v And b atmosphere (at least in 
the blue). This resembles our preliminary model atmosphere for 
HD 189733b ( jBerdyugTna 201 1\ with a high-altitude condensate 
layer which can reproduce the wavelength dependence of po- 
larization and albedo if average particle size is ~20nm. These 
clouds are highly reflective in the blue due to the small particle 



6. Geometrical albedos of hot Jupiters 

Detecting reflected light from exoplanets has been a challenging 
task since their discovery. First attempts were to extract a spec- 
tral signature from the observed spectrum of the star + planet 
system by comparing it to a modelled stellar spectrum, under 
the assumption that the reflected spectrum is a copy of the stellar 
spectrum (e.g., |Charbonneau et al.|[T9 99; Leigh et IT] |2003a] ). 
These provided upper limits on the optical reflected flux for 
some hot Jupiters. For transiting systems, detecting reflected 
light can constrain the geometrical albedo A g of the planet at a 
given wavelength. For instance, high-precision optical photom- 
etry in broad passbands with space telescopes MOST, CoRoT, 
and Kepler enabled detection or constrained secondary eclipses 
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Fig. 6. Geometrical albedo of hot Jupiters measured in various passbands (from Table [ij. Upper limits are shown by downward 
arrows. Some symbols are slightly moved from their central wavelength for clarity. In the upper panel, solid thick (blue) and dotted 
curves are the geometrical albedos of Neptune and Jupiter, respectively (Karkoschka, 1994), smoothed by a 100 nm boxcar. Dashed 
horizontal line indicates the albedo 0.791 of a semiinfinite conservative Rayleigh scattering atmosphere. Thin solid line is a modified 
Neptune albedo taking into account an additional absorption at 470-630 nm shown in the lower panel. 



and light curves effects from planetary phases. Modeling these 
light curves, while neglecting or evaluating the contribution due 
to planet's thermal radiation, resulted in estimates of albedos of 
transiting hot Jupiters. 

It is interesting to compare the polarimetric albedos of 
v And b and HD 189733b with those of other hot Jupiters. We 
collected known estimates of A„ in Table [3] together with other 
releveant data on the planets and host stars, as provided by the 
exoplanet.eu database and references therein ( |Schneider et al. 
201 l| l. The difference to our data is that passbands of the space- 
based measurements are significantly wider and centred in red 
wavelengths, and our albedo estimates limited to the range of 
350^-50 nm are unique. Also, measurements in several pass- 
bands for the same planet are still rare in the visible. 

This sample is rather homogeneous in terms of planetary and 
stellar parameters. In particular, all but one (HD189733b) plan- 
ets are orbiting late F - early G stars and at a distance smaller 
than 0.062 AU. Therefore, it makes sense to characterize this 
sample with average parameters: (T„) = 5930 K, (a) = 0.04 AU, 
(R p ) = lARj, (M p ) = 1.5Mj, (p) = 0.8 g sirT 3 , and <A g > = 0.2 
(average over all optical wavelengths). One can consider these 
averages as typical for a hot Jupiter in the vicinity of an early G 
dwarf. 

We find no significant correlations of the albedo with other 
planetary and stellar parameters. This can be in part due to large 



differences in the passbands. However, the wavelength depen- 
dence of the geometric albedo is of great interest, as it can well 
characterize the planetary atmosphere. 

In Fig. [6] we plot the available albedo estimates at the cen- 
tral wavelengths of the passpands. Upper limits are plotted as 
downward arrows. One can see a clear tendency for the albedo 
to decrease toward the red: from 0.6 in the near UV to essentially 
zero in the near IR. A few exceptionally high upper limits in the 
red are probably due to either very low signal-to-noise ratio or 
contamination by thermal emission. 

A comparison with the albedo of the Solar system giant plan- 
ets is educative. Earlier we concluded that the polarimetric albe- 
dos of the two hot Jupiters seem to be similar to the albedo of 
Neptune. The measurements by CoRoT and Kepler do not con- 
tradict this observation. Moreover, it is rather obvious that practi- 
cally all measurements strongly differ from those of the Jupiter's 
albedo, which is dimmer in the near UV and brighter in the red. 
Hence, Neptune's albedo is a more suitable reference for com- 
parison with hot Jupiters. 

As pointed out above, the albedo of the cold Neptune is 
mainly shaped by Rayleigh and Raman scattering on a high- 
altitude haze and H2 in the blue, and methan absorption in the red 
( |Sromovsky| [2005 ). In hot Jupiter atmospheres, H2 is still main 
scatterer, together with H 1, while absorption in the red is mainly 
due to TiO and H 2 (e.g., |Sudarsky et al.||2lX)0"l|Burrows etal) 
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Table 3. Geometrical albedos of hot Jupiters. 



Planet 




Passband 




M p 


P 


a 


T, 


R, 


Method + , 






[nm] 


IRA 


[Mj] 


[g cm -3 ] 


[AU] 


[K] 


[R B ] 


Reference 


HD 189733b 


0.62+0.31 


345-388 


1.178 


1.138 


0.86 


0.03142 


4980 


0.788 


Pol, [1] 




0.61+0.13 


390-490 


















0.28+0.16 


500-583 


















<0.26 


450-650 














Pol, [3] 


v And b 


0.53+0.27 


345-388 


1.35 


0.74 


0.37 


0.059 


6212 


1.631 


Pol, [2] 




0.67+0.24 


390-490 














" 




0.29+0.23 


500-490 
















HD209458b 


< 0.17 


400-700 


1.38 


0.714 


0.34 


0.04747 


6075 


1.146 


Ph, [4] 


CoRoT-lb 


0.20+0.08 


400-1000 


1.49 


1.03 


0.39 


0.0254 


5950 


1.11 


Ph, [5] 




<0.2 


560-860 














Ph, [6] 


CoRoT-2b 


0.16+0.03 


560-860 


1.465 


3.31 


1.31 


0.0281 


5625 


0.902 


Ph, [7] 


HATP-7b 


<0.18 


423-897 


1.421 


1.8 


0.78 


0.0379 


6350 


1.84 


Ph, [9] 




<0.13 


350-1000 














Ph, [8] 


TrES-2b 


<0.01 


423-897 


1.169 


1.253 


0.97 


0.03556 


5850 


1.0 


Ph, [10] 


TrES-3b 


<1.07 


550-700 


1.305 


1.91 


1.07 


0.0226 


5720 


0.813 


Ph, [11] 




<0.30 


700-850 


















<0.62 


850-1000 
















Kepler-5b 


0.12+0.04 


423-897 


1.431 


2.114 


0.90 


0.05064 


6297 


1.793 


Ph, [12,13] 


Kepler-6b 


0.11+0.04 


423-897 


1.323 


0.669 


0.36 


0.04567 


5647 


1.391 


Ph, [13] 


Kepler-7b 


0.32+0.03 


423-897 


1.614 


0.433 


0.13 


0.06246 


5933 


2.02 


Ph, [12,14] 


Kepler-8b 


<0.63 


423-897 


1.419 


0.60 


0.26 


0.0483 


6213 


1.486 


Ph, [12] 


Kepler- 17b 


0.10+0.02 


423-897 


1.31 


2.45 


1.35 


0.02591 


5630 


1.019 


Ph, [15] 


t Boo b 


<0.3 


466-499 


{1.2? 


>3.9 


2.80 


0.046 


6309 


1.331 


Sp, [16] 




<0.39 


385-611 


(1.2)* 












Sp, [17,20] 




<0.37 


590-920 


(1.2)* 












Pol, [21] 


HD75289A b 


<0.12 


400-900 


(1.6)* 


>0.42 


0.24 


0.046 


6120 


1.25 


Sp, [18] 




<0.46 


400-900 


(1.2)* 












Sp, [19] 



Notes. + Method to evaluate albedo: Pol - polarimetry, Ph - photometry, Sp - spectroscopy. * Assumed values. References: [1] Berdyugin a et al.| 
l |20Tl); [2] this paper; [3]|W iktorowic z1pOl39l > ; [4]|Rowe et aIT1([2008); [5]|Alonso et al.|p009l; [6]|Snellen et al.| ( |2009l; [7]|Snellen et al.|<|2010| >; 
[8]|Christiansen et al. 1(2010); [9]|Welsh et al.|(2010| >; [10]|Kipping & Spiegel|(201 1); [1 lJlWinn et al |(2008] >rU 2]|Kipping & Bakos)(201 It; [13] 
Desert et al. 1(201 la); [14~]|Demory et al.|(2011|>; [15]fD esert et al. H2011b|>; [16]|Charbonneau et al.|(199l?f ; [17] |Leigh et al |(2003afr : ?TTFf|Leigh 
et al.| ( |2003bV ; [19] |Rodler st al. 1 20081; [20] |Rodler et al.| ( |2010| i; [21] |Lucas et al.|q2009) >; 



2008 For tney et al.[[200^|Demory et al.|[20TT] l. These are per- 
haps the main opacities influencing the overall shape of the hot 
Jupiter albedo. Model predicts that the presence of high-altitude 
haze can significantly increase the albedo. This is revealed in 
HD 189733b with polarimetry (Berdyugina eTal] [20081 |20TT) , 
and transmission spectroscopy recently extended to the near UV 
confirms the overall presence of haze and Rayleigh scattering 



(|Pont et al.||2008[|Lecavelier d es Etang s et al.l|2008a||Sing et al. 



201 1) . This is in ine with our preliminary semi-empirical mode 
of HD 189733b (Berdyugina, 201 1 ) which explains the observed 
polarization by the presence of 20 nm dust condensates in its 
atmosphere. 

A closer inspection of the sample indicates that there is per- 
haps a systematic difference between exoplanets and Neptune, 
as shown in the lower panel of Fig. [6] It extends over about 
500-650 nm and reduces the albedo by ~0.2. This dimming 
can be caused by absorption in the sodium resonance doublet 
(590 nm) and the TiO electronic band systems (450-650 nm). 
These species play no role in the Neptune atmosphere but are 
essential for hot Jupiters. A varying strength of the absorption 
can strongly influence broad-band optical albedo values. For in- 



stance, the very low optical albedo of HD209458b ( |Rowe et al 
2008| > can be due to an extreme strength of this absorption fea- 



ture, which was also featured in the transmission spectroscopy 
by Sing et aL| ( |2008| >. They also detected an opacity increase to- 
wards shorter wavelenths which was interpreted as a signature of 
Rayleigh scattering (Lecavelier des Etangs et al. , 2008b). Hence, 



also for HD209458b, we can expect somewhat higher albedo in 
the blue. In other planets, such as Kepler-7b, a lack of sodium 
absorption and/or presence of haze could promote high albedo 
( |Demory et aL] |20TT) . 



7. Conclusions 

Our polarimetric study of v And has unveiled many unique and 
important implications for the dynamics and evolution of its 
planetary system and for physical properties of the planet b and 
hot Jupiters in general. 

- Determining the inclination of the u And b orbit and evalu- 
ating its radius, for the first time, allowed us to correct the 
planet mass, constrain its density, and reveal its inflated na- 
ture. 

- The found periastron longitude indicates that the apsidal res- 
onance known for planets c and d includes also planet b, 
which provides further constraints on their dynamical his- 
tory. 

- The limits obtained on the geometrical albedo idicate reflec- 
tive properties like in HD189733b. Both these hot Jupiters 
reflect most efficiently blue light and shine similar to 
Neptune in the Solar system. Their average optical albedo 
is close to that of Kepler-7b. It is possible that their high re- 
flectivity is due to high-altitude haze. 

- Combining the polarimetric albedos with those evaluated 
from optical spectroscopy and photometry for a sample of 
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planets, we obtained a unified picture of the geometrical 
albedo of an average hot Jupiter. Its albedo can be as high as 
0.6 in the near UV and close to zero at wavelengths longer 
than 700 nm. A broad absorption centred at 550 nm, prob- 
ably due to TiO and Nai, can strongly influence an average 
optical albedo of hot Jupiters. 

Our results define clear prospects for follow up studies of hot 
Jupiters. First, determinig albedo in the blue {A < 500 nm) for 
a larger sample of planets is crucial for clarifying whether the 
high reflectivity in the blue is a typical feature of hot Jupiters. 
Second, new self-consistent model atmospheres of hot Jupiters 
should be developed to explain such high reflectivity and to con- 
strain composition of their atmospheres, especially at high alti- 
tudes. Finally, polarimetric constraints on planet reflectivity are 
essential for separating reflected light from thermal emission in 
the flux measured through secondary eclipse photometry. 
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